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1 | INTRODUCTION 


Summary 

Equine coronaviruses (ECoV) are the only coronavirus known to infect horses. So 
far, data on ECoV infection in horses remain limited to the USA, France and Japan 
and its geographic distribution is not well understood. We carried out RT-PCR on 
306 nasal and 315 rectal swabs and tested 243 sera for antibodies to detect coron- 
avirus infections in apparently healthy horses in Saudi Arabia and Oman. We docu- 
ment evidence of infection with ECoV and HKU23 coronavirus by RT-PCR. There 
was no conclusive evidence of Middle East respiratory syndrome coronavirus infec- 
tion in horses. Serological data suggest that lineage A betacoronavirus infections are 
commonly infecting horses in Saudi Arabia and Oman but antibody cross-reactivities 
between these viruses do not permit us to use serological data alone to identify 


which coronaviruses are causing these infections. 
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tracts and cause a remarkably diverse spectrum of diseases in 


Coronaviruses infect a wide range of mammals and birds, exhibit a humans and livestock. They typically have a restricted host range, 


marked tropism for epithelial cells of the respiratory and enteric infecting only their natural host and closely related animal species 


but do have the capacity to cross the species barrier to infect new 


tThese authors contributed equally to this manuscript and are co-first authors. hosts. Equine coronavirus (ECoV) is the only coronavirus known to 
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infect or cause disease in horses (Balasuriya, 2013). It belongs to the 
genus Betacoronavirus lineage A, as does human coronavirus OC43, 
bovine coronavirus and porcine hemagglutinating encephalomyelitis 
virus. Equine coronavirus was first isolated from faeces of a diar- 
rhoeic foal in 1999 (ECoV-NC99) in North Carolina, USA (Guy, Bres- 
lin, Breuhaus, Vivrette, & Smith, 2000), and was initially believed to 
only affect foals. Since 2010, there have been several reports of 
ECoV-associated respiratory and enteric infections in adult horses in 
Japan, Europe and the United States, but its global distribution is still 
poorly defined (Kooijman, Mapes, & Pusterla, 2016; Miszczak et al., 
2014; Oue, Morita, Kondo, & Nemoto, 2013; Pusterla, Holzen- 
kaempfer, Mapes, & Kass, 2015). 

Middle East respiratory syndrome coronavirus (MERS-CoV) is a 
betacoronavirus within lineage C first identified in humans in 2012 
and continues to pose a threat to global health (WHO, 2016 May 
27). The evidence points to dromedary camels as a natural host 
and a major source for human MERS-CoV infections (Alexandersen, 
Kobinger, Soule, & Wernery, 2014; Chu et al., 2015; Gutierrez, 
Tejedor-Junco, Gonzalez, Lattwein, & Renneker, 2015; Memish 
et al., 2014; Perera et al., 2013; Reusken et al., 2014). There is no 
convincing evidence of MERS-CoV infections in other domestic 
livestock species so far (Adney et al., 2016; Alexandersen et al., 
2014; Hemida et al., 2013; Meyer et al., 2015; Perera et al., 2013). 
Sequence similarity comparisons of dipeptidyl peptidase-4 (DPP4 
[CD26]), the functional receptor for MERS-CoV, have revealed that 
equine DPP4 is phylogentically closely related to human DPP4 and 
the binding affinity of MERS-CoV spike S1 domain for equine 
DPP4 is similar to that of human and camel DPP4, raising the pos- 
sibility that MERS-CoV infections may occur among horses (Barlan 
et al., 2014). Serological studies of horses originating from Spain 
and the United Arab Emirates gave negative results for MERS-CoV 
(Alexandersen et al., 2014; Meyer et al., 2015), but it is still rele- 
vant to investigate for evidence of MERS-CoV in horses in areas 
endemic for MERS-CoV. 

We have tested for coronavirus infections in apparently 
healthy horses in Saudi Arabia and Oman by testing nasal and 
rectal swabs by RT-PCR assays for MERS-coronavirus and using a 
pan-coronavirus RT-PCR with potential to detect all coronaviruses, 
hitherto known and unknown. Following the detection of ECoV 
and HKU23 coronavirus (HKU23) by the pan-coronavirus RT-PCR, 
we tested the swab specimens with specific RT-PCR assays for 
ECoV and HKU23 and tested the sera with serological assays to 
detect ECoV, bovine CoV (BCoV) (closely related to HKU23), as 
well as MERS-CoV. 


2 | MATERIALS AND METHODS 


2.1 | Specimen collection 


Serum samples (n = 243), nasal swabs (n = 306) and rectal swabs 
(n = 315) were collected from apparently healthy adult (>1 years of 
age) horses from East and Central provinces of the Kingdom of 


Saudi Arabia and Oman. Many horses had all three samples 


collected. These specimens were collected during the period Febru- 
ary 2014 to January 2015. The numbers of specimens collected from 
different geographic regions of Saudi Arabia and Oman are indicated 
in Table 1. The animal management practices and exposures to other 
domestic animal species such as dromedaries are indicated. Farming 
activities include activities related to date palm plantations as well as 
other farming activities. Racing activities refer to horse racing within 
Saudi Arabia rather than participation in international races. 


2.2 | Cells and viruses 


Vero cells (ATCC CCL-81) were used to culture and propagate 
MERS-CoV and the HRT-18G cell line (ATCC, CRL-11663, American 
Type Culture Collection, Manassas, VA, USA) for culture of ECoV 
and BCoV. The culture medium for both the cell lines was DMEM 
(Cat:10569; Life Technologies, Carlsbad, CA, USA) supplemented 
with 10% foetal calf serum (Cat:16140; Life Technologies). 

The prototype viruses used were MERS-CoV (EMC/2012 
obtained from Dr R Fouchier, Erasmus MC, Rotterdam), the ECoV- 
NC99 strain (obtained from Dr. Udeni B.R. Balasuriya, Maxwell H. 
Gluck Equine Research Center, University of Kentucky, USA) (Guy 
et al., 2000; Zhang et al., 2007) and the bovine coronavirus (ATCC 
BRCV-OK-0514-2). 


2.3 | RT-PCR testing 


Viral RNA was extracted from nasal and rectal swabs using the 
QlAamp viral RNA minikit (Qiagen, Hilden, Germany) following the 
manufacturer's instructions. RNA extracts were tested for evidence 
of conserved coronavirus nucleic acid genetic sequences using previ- 
ously reported RT-PCR assays (Chu et al., 2014), RTqPCR assay for 
MERS-CoV upE gene (Corman et al., 2012), RTgPCR assay for ECoV 
(Miszczak et al., 2014), and a RTqPCR assay for HKU23 reported 
below. RNA was reverse transcribed in a 20 ul reaction mixture con- 
taining 1x first-strand buffer, 5 mM DTT, 0.5 mM deoxynucleotide 
triphosphates (dNTP), 2.5 ng/yl random hexamers and 200 units of 
SuperScript Ill (Life Technologies). 

The pan-coronavirus nested PCR targeted the RNA-dependent 
RNA polymerase (RdRp) gene of coronaviruses developed by us as 
previously described (Chu et al., 2014). Using the cDNA synthesized 
as described above, first round PCR was carried out using forward 
primer 5'-GGKTGGGAYTAYCCKAARTG-3’ (position 15,287 of ECoV 
strain NC99; GenBank accession number EF446615) and reverse pri- 
mer 5/-TGYTGTSWRCARAAYTCRTG-3' (position 15,869 of ECoV 
strain NC99) and Phusion High Fidelity PCR Master Mix Kit (Thermo 
Scientific). A 25 pl reaction mixture contained 12.5 ul of 2x phusion 
master mix, 1.5 ul (10 pmol) forward primer, 1.5 pl (10 pmol) reverse 
primer, 4.5 ul ddH2O and 5 ul of cDNA. The PCR cycler conditions 
for the amplification were 98°C for 2 min, 40 cycles of amplification 
(98°C for 15 s, 48°C for 15 s, 72°C for 30 s), then 72°C for 2 min. 
An aliquot of 0.5 ul of the first PCR product was then further ampli- 
fied in a second round PCR prepared as first around PCR but using a 
new set of primers (forward primer 5’-GGTTGGGACTATCC 
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TABLE 1 Geographic location of specimens collected and animal management practices 


Animal 
Provider management Contact with Nasal Rectal Virus RNA 
Region City code practices other animals Animals in contact with swabs swabs Sera _ detected, if any 
Saudi Arabia Riyadh 1 Breeding Occasionally Camels, sheep, goat, chickens 14 14 14. ECoV(n=5) 
(Central) 
2 Farming, Occasionally Camels, sheep, goat (0) (@) 33 
Breeding 
Saudi Arabia Dammam 3 Farming, Frequent Camels, sheep, goat, chickens 103 103 19 HKU23 
(East) Breeding 
4 Breeding Frequent Camels, sheep, goat 6 6 6 
Saudi Arabia Al-Hasa 5 Farming, Occasionally Camels, sheep, goat 39 39 39 
(East) Breeding 
6 Farming Frequently Camels, sheep, goat, chickens 27 27 0 
7 Farming Frequently Camels, sheep, goat, chickens, dogs 26 26 15 
8 Farming, Frequent Camels, sheep, goat, chickens 22 22 12 
Breeding 
Saudi Arabia Al-Qassem 9 Racing Rare Camels, sheep, goat 20 19 19 
(Central) 
Saudi Arabia Qateef 10 Farming Frequent Camels, sheep, goat, chickens 41 41 24  HKU23 
(East) 
11, Farming Frequent Camels, sheep, goat, chickens i@) (@) 13 
Saudi Arabia Aljobail il Farming Frequent Camels, sheep, goat, chickens (0) ie) 5 
(East) 
Oman 13 Racing, Frequent Camels, dogs 8 18 44 
breeding 
Total 306 S15 243 


TAAGTGTGA-3’ (position 15,287 of ECoV strain NC99), reverse pri- 
mer 5/-CCATCATCAGATAGAATCATCAT-3’ (position 15,721 of 
ECoV strain NC99) which amplify a final PCR product of 440 bp. 
The final PCR amplicons (440 bp) were gel purified using QlAquick 
gel purification kit (QIAGEN) and analysed by DNA sequencing using 
the primers used for the second round PCR. To amplify an extended 
RdRp gene sequence of both equine CoV and HKU23-like viruses, 
PCR was prepared as above with forward primer 5’-TACGCAG 
GATGGTAATGCTG-3’ (position 14,767 of ECoV strain NC99) and 
reverse primer 5’-TTGTGCGCATTCATTCGCAAG-3’ (position 15,404 
of ECoV strain NC99) and 5 ul of cDNA as template for the reac- 
tion. PCR products with expected size of 658 bp were purified for 
DNA sequencing using the forward and reverse primers. 

We developed a two-step real-time quantitative RT-PCR (RT- 
qPCR) assay for the detecting HKU23 targeting N gene of the virus. 
Virus RNA was reverse transcribed as detailed above. Real-time PCR 
was performed using Power SYBR Green Master Mix (Life Technolo- 
gies) in a 20 wl reaction containing 5 pmol forward primer 5’- 
CTGCCACGATGGTATTTTTAC-3’ (position 29,758 of HKU23 strain 
362F) and 5 _ pmol primer 5’/-GGTATTGACAT 
CAGCCTGGT-3’ (position 29,852 of HKU23 strain 362F) with 40 
cycles of amplification (95°C for 15 s, 62°C for 40 s) in a 20 ul vol- 


ume followed by melting curve analysis. Positive result was inter- 


reverse 


preted by a detectable Ct value >40 considered together with a 
amplicon melting temperature within the range of 81.64°C to 
82.64°C. 


The details of all RT-PCR and real-time PCR assays used in this 


study are summarized in Table S1. 


2.4 | PCR amplicon sequencing and analysis 


The primers used for sequencing were those used to generate the 
PCR products. DNA amplicons were purified and sequenced in both 
forward and reverse directions with the PCR primers using BigDye 
Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Foster 
City, CA) and a 96-capillary 3730xl DNA Analyzer (Applied Biosys- 
tems) according to manufacturer's instructions. Sequences were 
assembled using BioEdit 7.0 and MUSCLE for multiple sequence 
alignment. Phylogeny analysis was performed by maximum likelihood 
method using PhyML with GTR nucleotide substitution model and 
SH-like approximate likelihood ratio test for estimating branch 


support. 


2.5 | Serology 


All sera were tested using previously validated MERS-CoV spike pseu- 
doparticle neutralization test (ppNT) (Hemida et al., 2014; Perera 
et al., 2013), and the positive sera were confirmed using microneutral- 
ization (MN) test and plaque reduction neutralization tests (PRNT) 
using live MERS-CoV strain EMC in a Biosafety Level 3 containment 
laboratory (Park et al., 2015). Selected sera were tested for antibodies 
to ECoV and BCoVs (Perera et al., 2013). The North American ECoV- 
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NC99 strain was propagated in HRT-18G cells to make virus stocks 
for MN assay. ECoV MN assays were carried out by mixing 100 ul of 
serial 2-fold dilutions (1:10 to 1:640) of heat-inactivated sera (56°C for 
30 min) with equal volumes of 200 tissue culture infectious doses 
(TCIDs0) of ECoV in 96-well microtitre plates (Cat:?2096, TPP tissue 
culture plates, Switzerland). Plates were incubated for 1 hr at 37°C. 
The virus-serum mixture was then added in quadruplicate to HRT- 
18G cell monolayers in 96-well microtitre plates followed by 1 hr 
adsorption at 37°C. The plates were incubated at 37°C in 5% COz ina 
humidified incubator for 5 days, when virus back titrations were 
inspected to confirm that the virus dose used was within the expected 
range. The antibody titre was the highest serum dilution that resulted 
in complete neutralization of cytopathic effect in >2 of the four wells 
replicate wells at each serum dilution. Antibody titre refers to the dilu- 
tion of the serum in the serum dilution—virus mixture. All antibody 
titres are denoted as the reciprocal of the titre. 

BCoV is ubiquitous in cattle worldwide and is genetically and 
antigenically related with HKU23, a coronavirus endemic in dromed- 
aries (Woo et al., 2014, 2016). Antibody to BCoV was tested using 
MN assay as described previously (Park et al., 2015; Perera et al., 
2013). Serology for BCoV and ECoV was carried out with culture 
medium without foetal calf serum. Known positive and negative con- 


trol sera were included in all serology assays. 


2.6 | Reference sera used in investigation of cross- 
reactivity in serological assays 


Cross-reactive responses in the MERS-CoV ppNT, BCoV MN and 
ECoV MN assays were assessed using a panel of previously charac- 
terized immune or hyperimmune reference antisera generated 
against alpha, beta and gamma coronaviruses. These sera were either 
obtained from BEI Resources (animal CoV reagents supplied to BEI 
by Dr Linda Saif http://www.beiresources.org/About/BEIResources.a 
spx) or kindly provided by Dr. Linda Saif or Dr Stanley Perlman. 
Equine sera known to be seropositive to ECoV were provided by Dr. 
Nicola Pusterla, Department of Veterinary Medicine and Epidemiol- 
ogy, School of Veterinary Medicine, University of California, Davis, 


CA, USA. The homologous antibody titres to the immunized virus 


TABLE 2 RT-qPCR-positive specimens 


Pan CoV nested 
RT-PCR result 


were also obtained, where available, from the providers of these 
respective antisera. The serum panel for alpha coronavirus included 
sera generated against porcine respiratory coronavirus, feline infec- 
tious peritonitis virus, canine enteric coronavirus and porcine trans- 
missible gastroenteritis virus. Sera produced against mouse hepatitis 
virus strains JHM and A59, SARS coronavirus, BCoV, ECoV repre- 
sented the betacoronavirus group. Furthermore, an immune serum 
for infectious bronchitis virus was used as representative for the 
gammacoronavirus group. All these reference sera, as well as addi- 
tional ECoV immune sera and MERS-CoV-negative camel sera, were 
tested in ppNT or MN assays in parallel against MERS-CoV, BCoV, 
as well as ECoV, to assess possible cross-reactive antibody 
responses among these three coronaviruses. Some of these refer- 
ence sera, with the exception of horse and camel sera, had been 
previously tested against MERS-CoV and BCoV, and the results were 
reported by us (Hemida et al., 2014). 


3 | RESULTS AND DISCUSSION 


All nasal and rectal swabs were tested using the pan-coronavirus 
nested RT-PCR assay, ECoV and MERS-CoV upE gene RT-qPCR 
assays. All swabs were negative in the MERS-CoV RTqPCR assay. 
Five rectal swabs (5/315) were positive for coronavirus RNA by the 
pan-coronavirus assay (Tables 1 and 2). Sequencing of the 440 base 
pair nucleotide fragment from the RdRp gene confirmed four of 
these to be ECoV and the other to be dromedary camel HKU23-like 
virus designated as equine HKU23 CoV SA1354. Additional 
sequence was obtained from these specimens to generate a nucleo- 
tide sequence of 933 nucleotides (position 14,780 to 15,712 of 
ECoV strain NC99) which was used for analysis of nucleotide differ- 
ences in comparison with other coronaviruses (GenBank accession 
numbers KY458228-KY458232) (Table 3) and phylogenetic analysis 
which confirmed the identification of the coronaviruses to be ECoV 
and HKU23, respectively (Figure 1). This result led to the testing of 
306 nasal and 315 rectal swab specimens by specific RTqPCR assays 
for ECoV and RTgPCR assay for HKU23. Testing by the ECoV speci- 


fic RTqPCR led to the detection of five rectal swabs specimens, four 


Location and with sequencing MERS-CoV EqCoV_RT- HKU23 RT- HKU23 RT- 

Swab type provider code of PCR product UpE_RT-qPCR qPCR_Ct qPCR Ct value qPCR Tm? 
Nasal Al-Quassem (10) Negative Negative Negative 37.46 82.26 
Rectal Dammam (3) HKU23 Negative Negative 34.67 82.25 
Rectal Riyadh (1) Equine CoV Negative 35.99 Negative 

Rectal Riyadh (1) Equine CoV Negative 29.38 Negative 

Rectal Riyadh (1) Equine CoV Negative 30.63 Negative 

Rectal Riyadh (1) Negative Negative 36.32 Negative 

Rectal Riyadh (1) Equine CoV Negative 34.87 Negative 


*Tm of camel HKU23 is 82.138 degree C. Screening horse samples for HKU23 RTqPCR allowed Tm range 81.638 to 82.638 degree C. 
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Equine CoV SA1360 
Equine CoV SA1362 
Equine CoV SA1364 
Equine CoV SA1382 


B 


HCoV 0€43 (AY391777) 


Bovine CoV strain Mebus (U00735) 


Equine HKU23 CoV SA1354 


Bovine CoV DB2 1983 (DQ811784) 


0.09 


Equine CoV strain NC99 (EF446615) 

Equine CoV RNA strain Tokachi09 (LC061272) 

Equine CoV RNA strain Obihiro12-2 (LC061274) 

Equine CoV RNA strain Obihiro12-1 (LC061273) 

[Porcine hemagglutinating encephalomyelitis virus strain VWW572 (DQ011855) 


Dromedary camel CoV HKU23 strain HKU23-368F 21-May-13 UAE Dubai (KF906251) 
Bovine CoV E-DB2-TC 30-Nov-96 USA Ohio (FJ938063) 


HCoV HKU1 (NC 006577) 


Equine coronavirus 


HCoV 0C43 serotype 0€43-Paris (AY585229) 
HCoV 0C43 strain ATCC VR-759 (AY585228) 
HCoV 0€43/human/USA/851-15/1985 8-Jan-85 USA Nashville TN (KF530060) 
HCoV 0€43/human/USA/873-16/1987 12-Mar-87 USA Nashville TN (KF530077) 
HCoV 0€43/human/USA/873-19/1987 17-Mar-87 USA Nashville TN (KF530083) 
HCoV 0€43/human/USA/872-5/1987 10-Feb-87 USA Nashville TN (KF530086) 
Canine respiratory CoV strain K37 2008 lung South Korea (JX860640) 
Bovine CoV genomic RNA strain Kakegawa (AB354579) 
Bovine CoV strain Quebec (AF220295) 


Human coronavirus OC43 


Human enteric CoV 4408 1988 Germany (FJ415324) 

Human enteric CoV strain 4408 17-Apr-09 (FJ938067) 

BetaCoV 1 isolate Buffalo CoV B1-24F 1-Oct-14 Bangladesh (KU558922) 

BetaCoV 1 isolate Buffalo CoV B1-28F 1-Oct-14 Bangladesh (KU558923) 
Dromedary camel CoV HKU23 strain HKU23-265F 3-Apr-13 UAE Dubai (KF906249) 


Bovine coronavirus/HKU23 cluster 


Bovine CoV strain BCV-AKS-01 30-May-15 China Xinjiang (KU886219) 
White-tailed deer CoV US/OH-WD470/1994 1-Jan-94 USA Ohio (FJ425187) 

Bovine respiratory CoV bovine/US/OH-440-TC/1996 1-Jan-96 USA Ohio (FJ938066) 
Sambar deer CoV US/OH-WD388/1994 1-Jan-94 USA Ohio (F1425189) 

Sambar deer CoV US/OH-WD388-TC/1994 1-Jan-94 USA Ohio (FJ425188) 
Waterbuck CoV US/OH-WD358-GnC/1994 1-Jan-94 USA Ohio (FJ425185) 

Waterbuck CoV US/OH-WD358-TC/1994 1-Jan-94 USA Ohio (FJ425184) 

Waterbuck CoV US/OH-WD358/1994 1-Jan-94 USA Ohio (FJ425186) 


FIGURE 1 Phylogenetic analysis of equine coronavirus, human coronavirus OC43 and bovine coronavirus/dromedary camel coronavirus 
HKU23 cluster. The maximum likelihood phylogenetic tree rooted to HKU1 as an outgroup was generated using virus RdRp gene alignment of 
934 nt by PhyML with approximate likelihood ratio test (aLRT) SH-like branch support. Sequences generated from this study are highlighted in 


red. GenBank accession numbers of are indicated in brackets. 


TABLE 4 MERS-CoV serology results 


Uenone Reciprocal antibody titre 


Horse ID provider code MERS-ppNT MERS MN 
5-36 Al-Hasa (5) 40 <20 
5-41 Al-Hasa (5) 40 <20 
9-2 Al-Qassem (9) 160 20 
13-12 Oman (13) 40 <20 
2-12 Riyadh (2) 40 <20 


MERS PRNT90 MERS PRNT50 BCoV MN ECoV MN 
<20 20 <20 80 
<20 20 <20 20 
40 80 <20 <20 
<20 20 <20 <20 
<20 20 20 40 


ppNT, pseudotype neutralization test; MN, microneutralization test; PRNT, plaque reduction neutralization test. 


of which had been detected by the pan-coronavirus RT-PCR assay 
as well as an additional positive rectal swab (Ct value 36.3) from a 
horse in the Riyadh area (Table 2). Using the HKU23-specific 
RTqPCR, two positives were detected including SA1354 identified as 
HKU23 by sequencing the pan-coronavirus RT-PCR assay amplicon 
and an additional positive result with a high Ct value of 37.6 from a 
nasal swab sample. None of the other nasal swabs had evidence of 
coronavirus RNA. Compared to the pan-coronavirus assay, the speci- 
fic ECoV- and HKU23-specific assays detected one additional speci- 
men that was weakly positive for each of these viruses (Table 2). 
The five pan-coronavirus PCR CoV-positive specimens had Ct 
ranging from 29.4 to 35.9 in either the ECoV or HKU23 specific 
assays. The pan-coronavirus RT-PCR was not sensitive at detecting 
specimens with Ct higher than 36 in the ECoV or HKU23 assays. 


However, the results do demonstrate the importance of using the 
pan-coronavirus assay when investigating coronavirus activity in ani- 
mal species in new geographic areas. Our initial focus was to explore 
for evidence of MERS-CoV infections in horses in Saudi Arabia and 
Oman. None of the rectal swabs were positive for MERS-CoV RNA, 
but the application of the pan-coronavirus RT-PCR to these speci- 
mens revealed the presence of two other coronaviruses infecting 
horses in this region, that is ECoV and HKU23. 

All seven horses from whom ECoV or HKU23 were detected 
were apparently healthy. The specimens positive for ECoV RNA all 
came from different horses at a horse breeding facility in Riyadh. 
The two horses positive for HKU23 RNA were from Dammam and 
Qateef, Saudi Arabia, where horses were used for farming, and in 


Dammam, also for breeding (Tables 1 and 2). In both cases, these 
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TABLE 5 Cross-neutralization titres (denoted as reciprocal titres) for Middle East respiratory coronavirus (MERS-CoV), bovine coronavirus 
(BCoV) and equine coronavirus (ECoV) in hyperimmune or naturally infected sera known to be positive for different coronaviruses 


Sera panel—Catalogue numbers provided 


Genus for sera obtained from BEI Resource 


Alpha coronavirus NR2727 — Guina pig antiserum to canine 


coronavirus 


NR460 — pig antiserum to porcine respiratory 
coronavirus 


NR2518 — Guina pig antisera feline infectious 
peritonitis virus 


NR458 — pig antisera for porcine transmissible 
gastroenteritis virus 


Beta coronavirus ECoV-positive sera 459 (horse, USA) 
ECoV-positive sera 463 (horse, USA) 
ECoV-positive sera 472 (horse, USA) 
ECoV-positive sera 486 (horse, USA) 
ECoV-negative sera 2 (horse, Hong Kong) 
BCoV-positive camel serum #740293° 
BCoV-positive camel serum #467468° 
C99-10 — BCoV antisera from guinea pig 


Bcov antisera from germ free bovine calf 


NR456 — BCoV antisera from gnotobiotic calf 


NRC772 — Rabbit antisera SARS S protein high titre 


Mouse hepatitis virusJHM strain) hyperimmunized 
mouse dam 1 


Mouse hepatitis virusJHM strain) hyperimmunized 
mouse dam 2 


Mouse hepatitis virus(A59 strain)-infected mice 
NRC774 — antisera for SARS coronavirus sero titre 


NRC769 — Rabbit antisera to SARS S protein sero 
titre 


Gamma coronavirus NR2515-Guina pig antiserum to infectious 


bronchitis virus 


Homologous antibody titre by MERS- MERS BCoV’ ECoV 
ELISA unless specified ppNT MN MN MN 
4,094? <20 <20 <20 <20 
1,200° <20 <20 <20 <20 
2,000° <20 <20 <20 <20 
1,400? <20 <20 <20 <20 
qPCR negative; ELISA OD 1.4597 <20 <20 160 >1,280 
aPCR positive; ELISA OD 1.5257 = <20 <20 1,280 >1,280 
qPCR negative; ELISA OD 1.458° 20 <20 640 >1,280 
qPCR negative; ELISA OD 1.939° 20 <20 1,280 >1,280 
Not relevant <20 <20 <20 <20 
Not relevant <20 ND 640 80 
Not relevant <20 ND 80 <20 
20,480° <20 <20 320 <20 
580° — neutralisation titre <20 <20 1,280 80 
10,000° <20 <20 80 <20 
640 <20 <20 <20 <20 
1,778° — neutralisation titre <20 <20 <20 <20 
363° — neutralisation titre <20 <20 <20 <20 
1,000‘ — neutralization titre <20 <20 <20 <20 
<10 <20 <20 <20 <20 
<10 <20 <20 <20 <20 
50,000? <20 <20 <20 <20 


ppNT, Pseudoparticle neutralization test; MN, micro neutralization test; PRNT, plaque reduction neutralization test; BCoV, bovine coronavirus; ECoV, 


equine coronavirus; OD, optical density. 


Equine antisera against ECoV were kindly provided by Professor Nicola Pusterla, Department of Veterinary Medicine and Epidemiology, School of 


Veterinary Medicine, University of California, Davis, CA, USA. 


Homologous titres or OD results obtained from BEI resources, Linda Saif, Stanly Perlman and Udeni Balasuriya represented as 


ab & d respectively. 


“Field collected dromedary sera from Kazakhstan, a region known to be free of MERS-CoV infection. 


horses frequently came into contact with camels and other domestic 
livestock. HKU23 coronavirus is known to be endemic in camels 
(Woo et al., 2014). 

We investigated for serologic evidence of MERS-CoV infections 
in the 243 horse sera using ppNT screening followed by endpoint 
titrations. Any sera positive at a screening reciprocal antibody titre 
of 20 in ppNT were confirmed using MN assay and PRNT assay 
using live MERS-CoV strain EMC. Possible serologic cross-reactions 
with ECoV and BCoV were evaluated in MN assays which are the 
most specific serological assays available. Of the 243 sera tested, 
only five (2%) had detectable MERS-CoV antibody in the MERS- 
ppNT assay at a screening titre of 20; four of these were from Saudi 
Arabia, and one was from Oman. The endpoint MERS-ppNT titres of 


these five sera ranged from 40 to 160 but only serum 9-2 which 


had the highest MERS-ppNT titre was positive in the MERS-CoV 
MN test (titre 20) (Table 4). In the confirmatory PRNT which is con- 
sidered the “gold standard” serological test for MERS-CoV, horse 
serum 9-2 reduced plaque numbers by 90% (PRNTo9o) up to a titre 
of 40. This serum was negative in ECoV and BCoV MN assays 
(Table 4). This serum was from a racing horse stable in Al-Qassem 
and may represent a rare example of transmission of MERS-CoV to 
horses or cross-reaction with another yet undocumented coronavirus 
infecting horses. This horse was not known to frequently come into 
contact with camels. Three of the other sera (serum ID nos 5.35; 
5.41; 2.12) positive in the MERS-CoV ppNT assay had detectable 
ECoV MN antibody and one of these also had detectable BCoV MN 
antibody (Table 4). Thus, these may well represent evidence of 


cross-reactive responses. 
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TABLE 6 Seropositivity to BCoV and ECoV in horses from different geographic locations 


Provider code and Reciprocal Reciprocal 
Location management system Sample ID BCoV MN titre ECoV MN Titre PCR results 
Riyadh 1 (Breeding) 1-1? <20 20 ECoV PCR positive 
1-2? <20 20 ECoV PCR positive 
v1-3° 20 160 ECoV PCR positive 
1-4 <20 80 
1-5 <20 160 
1-6 <20 40 
1-8? <20 640 ECoV PCR positive 
1-12° 80 320 ECoV PCR positive 
2 (Farming, breeding) 2-13 40 160 
2-14 <20 80 
Dammam 3 (Farming, breeding) 3-11 <20 <20 
3-12 <20 40 
3-41? <20 <20 HKU23 PCR positive 
Al-Hasa 5 (Breeding) 5-56 <20 40 
9-5/ 160 320 
5-58 40 80 
5-59 80 160 
8 (Farming, breeding) 8-2 20 80 
8-3 1:80 40 
8-4 <20 640 
8-5 <20 <20 
8-6 <20 20 
7 (Farming) 7-1 <20 20 
7-2 20 <20 
7-3 <20 <20 
7-4 <20 20 
7-1 <20 <20 
7-2 <20 <20 
Al-Qassem 9 (Racing) 9-A3 20 <20 
9-A4 <20 40 
9-A5 <20 <20 
9-B1 <20 <20 
o-E1l 20 20 
9-C2 80 640 
Geel <20 <20 
9-E2 <20 <20 
Qateef 10 (Farming) 10-13 <20 80 
10-14 <20 20 
10-15 <20 80 
10-16 20 160 
11 (Farming) 11-8 80 40 
11-9 <20 320 
11-11 <20 <20 
11-12 <20 <20 
11-13 <20 80 


(Continues) 
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TABLE 6 (Continued) 
Provider code and Reciprocal Reciprocal 
Location management system Sample ID BCoV MN titre ECoV MN Titre PCR results 
Aljobail 12 (Farming) 12-1 40 40 
Oman 13 (Racing, breeding) 13-1 <20 20 
= 2 <20 80 
13-3 <20 40 
13-4 20 160 
13-5 20 40 
13-6 20 160 
13-10 20 40 
i= Ala <20 80 


*Horses found to be coronavirus RNA positive in swabs by RT-PCR or RT-qPCR. 


The extent of the serological cross-reactivity of MERS-CoV, 
ECoV and BCoV was systematically investigated. We did not have 
an isolate of HKU23 virus available for serological testing and BCoV 
was used as a virus that is genetically closely related to HKU23 with 
pairwise amino acid similarity of 94.1% in the spike protein gene, 
the determinant of antigenic cross-reaction in neutralization tests 
(Woo et al., 2014). Horse sera obtained from the USA reported to 
be strongly reactive by ELISA (optical density ranging from 1.458 to 
1.939) for ECoV showed high titre (>1,280) for ECoV in MN assay 
as expected, but also had high antibody titres (ranging from 160 to 
1,280) in BCoV MN tests (Table 5). As these horse sera were from 
horses in the “field,” they may well have been exposed to multiple 
coronaviruses. Thus, these results may reflect cross-reactivity or co- 
infection of US horses with ECoV- and BCoV-like viruses. Two of 
these sera showed low titres (20) in the MERS-CoV ppNT assay but 
were negative in the more stringent MERS-CoV MN assay. As 
MERS-CoV is not circulating in the USA, these results suggest cross- 
reactivity between MERS-CoV and ECoV or related CoVs may occur, 
albeit at low titre. Similarly, a BCoV immune serum from an experi- 
mentally infected gnotobiotic calf showed detectable, but 16-fold 
reduced antibody titre with ECoV but no cross-reaction with MERS- 
CoV. Spike proteins of equine coronavirus and bovine coronavirus 
have an overall amino acid similarity of 80.8% and 73.1% similarity 
of the $1 region which contains the receptor-binding domain and 
neutralization epitopes. This probably explains the large difference in 
neutralization titres observed in Table 5. A dromedary serum from 
Kazakhstan, a region known to be free of MERS-CoV activity 
(Miguel et al., 2016), had detectable antibody to BCoV (probably 
reflecting HKU23 or BCoV infection) and also to ECoV. We cannot 
conclude whether this reflects cross-reactivity or infection with 
BCoV, HKU23 and/or ECoV. None of the beta coronavirus immune 
sera to SARS-CoV or mouse hepatitis virus, nor immune sera to 
alpha or gamma coronaviruses gave cross-reactions in BCoV or 
ECoV MN assays (Table 5). As had been previously reported, none 
of these dromedary sera cross-reacted with MERS-CoV (Hemida 
et al., 2014). 

Fifty-four horse sera from Saudi Arabia and Oman were selected 


to represent different collection locations, excluding those five sera 


previously reported to be MERS-CoV ppNT positive (Table 6). These 
sera were seronegative for MERS-CoV (as expected); 40 (74%) of 
them had detectable MN antibody to ECoV (titres ranging from <20 
to 640; median 40), and 18 (33.3%) had MN antibodies to BCoV 
(titres ranging from <20 to 160; median <20). The reciprocal geomet- 
ric mean antibody titres (assigning sera with titres <20 the value of 
10 for computational purposes) for ECoV and BCoV were 44 and 
25.5, respectively. There were twenty-four sera with undetectable 
antibody to BCoV with detectable ECoV antibody titres ranging from 
20 to 640. The titres to ECoV were higher than or equal to titres to 
BCoV, with two exceptions where the titre to BCoV (20) was mar- 
ginally higher than that with ECoV (Table 6). These results are com- 
patible with the RT-PCR detection of ECoV in rectal swabs in horses 
and probably suggest that ECoV or a closely related virus is com- 
monly infecting horses in many regions of the Arabian Peninsula. 
Two of the five horses in whom ECoV RNA was detected by 
RTqgPCR had low ECoV MN antibody titres (20) while the other 
three had titres ranging from 160 to 640 (Table 6). This may reflect 
virus shedding persists into convalescence or that re-infection of 
previously seropositive horses can occur. 

The RTqPCR detection of two HKU23-like viruses indicates that 
HKU23 which is common in dromedaries in the region is also infect- 
ing horses. In the absence of an HKU23 isolate for use in neutraliza- 
tion tests and given the potential serological cross-reactivity 
between these three viruses, it is difficult to conclude how com- 
monly equines are infected with HKU23. The antibody titres to 
BCoV, which is genetically closely related to HKU23 in the spike 
protein gene, may suggest that BCoV- or HKU23-like virus infection 
of horses does occur in Saudi Arabia. 

The serological evidence of one of the 243 sera with detectable 
MERS-CoV antibody titres in the PRNT99 test raises the question 
whether MERS-CoV may occasionally cross species to horses. This 
serum did not have any detectable antibody to ECoV or BCoV. 
However, it is not possible to exclude the possibility that this obser- 
vation reflects cross-reaction with ECoV, BCoV or other yet 
unknown coronaviruses. Confirmation of natural infection of horses 
with MERS-CoV can only come from detection of MERS-CoV RNA 


in horse specimens. Recent experimental infection of horses with 
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MERS-CoV did not result in virus replication or seroconversion 
(Adney et al., 2016). 

In conclusion, the serological data suggest that lineage A 
betacoronaviruses are commonly infecting horses in Saudi Arabia 
and Oman but antibody cross-reactions between these viruses do 
not permit us to use serological data alone to identify which 
coronaviruses are causing these infections. RT-PCR detection of 
ECoV and HKU23 in equine swabs confirms the circulation of 
these two viruses in horses in Saudi Arabia. Cocirculation of 
related viruses may provide potential for recombination, a poten- 
tial means of generating genetic diversity and facilitating host 
jumps in coronaviruses. This is the first report of HKU23 being 
detected in horses and the first detection of ECoV in Asia outside 


of Japan. 
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